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We study population dynamics of a trapped spin-1 Bose gas above the Bose-Einstein transition 
temperature. Starting from the semiclassical kinetic equation for a spin-1 gas, we derive a coupled 
rate equations for the populations of internal states. Solving the rate equations, we discuss the 
dynamical evolution of spin populations. We also estimate the characteristic timescale in which the 
system reaches equilibrium. Finally, we briefly discuss how the presence of the condensate will affect 
the population dynamics. 

I. INTRODUCTION 

Bose-Einstein condensation (BEC) in dilute atomic gases have been extensively studied both theoretically and 
experimentally since its experimental realization in 1995 In particular, BEC involving internal degree of freedom 

has been extensively studied. The first experiment confining BEC with multi internal states was achieved by JILA 
group 3. Shortly afterward, they succeeded in confining ^^Rb atoms with two internal states \F = 2,mp = 1) and 
\F = 1,TOF = —1), which is called spin- 1/2 Bose gas Spin- 1/2 Bose gases are known to exhibit the collective spin 
dynamics due to the exchange effect even above the Bose-Einstein transition temperature (Tbec) 0, & 0, B S [13, EI ■ 

On the other hand, the spinor condensate trapped in an optical trap has been achieved experimentally. An optical 
trap can trap atoms with maintaining spin degree of freedom in contrast to magnetic trap where the spin degree of 
freedom is frozen. The MIT group first succeeded in creating a spinor condensate in an optical trap for ^'^Na (so 
called spin-1 BEC) [T^- In the spin-1 system, the spin-spin interaction that describes the collision of two atoms 
exchanging there internal states play s an important role, even though they are very small compared with the usual 
spin-independent interaction [H, 111 [13, [H, M, [13, [ll, [ll, HO, [H ■ 

Moreover in the spin-1 system, these components couple automatically through intrinsic spin-spin interaction. This 
is in contrast to the spin- 1/2 system, where these components couple with each other through the exchange effect. 
We thus expect the spin-1 Bose gas to exhibit much richer spin dynamics than in the spin-1/2 system. 

Even though most of the studies in the dynamics of spin-1 Bose gases have focused on a pure condensate, there 
have also been a growing interest in finite-temperature properties fisl . [l9l [2^ . In our previous paper, we discussed 
the spin- wave collective modes with dipole symmetry above Tbec ^M- We showed that the spin-1 Bose gas exhibits 
a fascinating dynamics even if above Tbec a-nd we indicated that the spin-1 Bose condensate gas below Tbec will 
exhibit much richer dynamics due to the interaction between condensate and noncondensate atoms in addition to the 
spin-spin interaction. 

One of the most interesting experiments of the dynamics of spin-1 Bose condensed gases at finite temperatures is the 
population dynamics of each internal state [l9| . In this experiment, one found that the condensates do not conserve 
the total spin by themselves. In addition, the existence of the condensate enhances equilibration of the noncondensate 
component [2^. We are interested in this property. 

In this paper, we first derive rate equations for spin-1 Bose gases above Tbec- Solving the rate equations, we 
discuss the population dynamics. We estimate characteristic time required to reach equilibrium. In the experiment, 
the time evolution of the particle numbers of each internal state was observed below Tbec where both condensate 
and noncondensate atoms are present [l^ . Comparing our theory, which only includes noncondensate atoms, with 
this experiment, we discuss how the existence of the condensate atoms will affect to the population dynamics. 
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II. SPIN-1 KINETIC EQUATION 



We consider a gas of bosonic atoms with spin F = 1 in an optical trap, so called spin-1 Bose gas. Each atom has 
three hyper-fine spin states; mp= 1,0, —1. In second quantized form, the Hamiltonian for this system is given by 

H = J2 f dr^ (r,t) I ^(r) | j)^, (r, i) 

+ f E / (r, t) *] (r, t) (r, t) ^, (r, t) 



where ^'J (r, t) is the Bose field operator of the Heisenberg expression satisfying the equal-time commutation relation: 
(r, , ^1 (r', t) = 6ijS {r~r'), and i,j indicate the hyperfine spin states mp. The hat indicates a second 
quantized operator. Here. Hq (r, t) is the single atom Hamiltonian given by 



Hoir,t) = 



2m 



-gfiB (r, t) r + E £ 



a(3 



(2) 



Q/3 



where the double underline indicates a 3 x 3 matrix, and a,/3 indicate spin components; x,y,z. The first term of 
Eq.(l2]) describes the center-of-mass motion of an atom in a harmonic trap V {r) = ^ iyt^'^ + '^'yV^ + '^'z^'^')i where 
m is the mass of a single atom. The second term describes the linear Zeeman energy, where B°' is spin-a component 
of the magnetic field, and 5" {a = x,y,z) denotes the spin-1 matrix of spin-a component, and g is Lande g factor 

and fiB is Bohr magneton. The third term is the quadratic Zeeman energy, in which we introduced the quadratic 

magnetic field B^l^ and the quadrupole operator Q"'^ is given by Q"^^ = (l - ^3^^) + - Sap'j f The 

coupling constants go and 52 are given by go = (47r?i^/3m) (qq + 202) and g2 ~ (47r7i^/3m) (02 — ao), qq and 02 being 
the s-wave scattering lengths for collision channels with total spins T = Q and 2, respectively. In this system, only 
coupling process between states with = and = ±1 can be permitted. 

The time evolution of the nonequilibrium system is described by the Wigner distribution function 



M^y(p,r) = y"drVP-'-'/''(*](r + r72)^',(r-r72)), 



(3) 



Knowledge of this function allows one to calculate various nonequilibrium physical quantities, such as the local density 
given by 



^(r, t) 



dp 



(27r?i)' 



:^(p,r,t) 



(4) 



where rnj (r, = (r, t) (r, ^))- Using a semiclassical approximation to describe atomic motion in terms of a 
phase-space distribution function, we obtain the semiclassical kinetic equation p^ : 



h 



dt m 

where the effective potential is described as 



(5) 



(6) 



We defined the magnetization as (r, t) = Tr ^(r, t) . The collision integral on the right hand side of Eq.([5|) 
describes collisions between atoms, whose explicit expression is given in Ref.(2^ 



3 



Throughout this paper we assume relatively high temperature regime where the phase space density is sufficiently 
low, i.e. Wij <^ 1. Under this assumption, in the collision integral 1^, we only retain terms to second order in Wij. 
This is equivalent to neglecting quantum degeneracy in the center of mass motion of atoms. 

The static thermal equilibrium distribution is determined from the condition = 0. This leads to the Maxwell- 
Boltzmann distribution 

W°=S.,W^ = e-''^(^^-'^'\ (7) 
where is the single-particle excitation energy 

ti = -^ + Uu, (8) 
and is the chemical potential of the i state. The condition for the equilibrium chemical potentials is given by 

Hi + ^-1 = 2^0- (9) 



III. RATE EQUATIONS FOR THE POPULATION DYNAMICS ABOVE THE BOSE-EINSTEIN 

TRANSITION TEMPERATURE 

We consider the time evolution of atoms in each internal state above the transition temperature. The particle 
number of the internal state i is given by 

iV,; {t) = j dv j ^r^Wu (P, r, t) . (10) 

The equation for the particle number is obtained by integrating kinetic equation Eq.([5|) over r and p. In this section, 
we assume that the system is close to equilibrium at each time, and thus approximate the distribution function with 
the Maxwell-Boltzmann form: 

W^j (p, r, t) = S,jf, (p, r, t) = % exp {- [e, (p, r, t) - fi, (t)] /ksT, (t)} , (11) 

where Hi and depend on time. With this approximation, the collision integral reduces to 



x2gl {s°^,Sil.^Sf^,S^,^ + S'^,S^,^S^,Sil.^ [U (Pa) ff (pd - h (p) (P2)] ■ (12) 



We note that the terms proportional to and .90.92 hi the collision integral vanish because of conservation of momen- 
tum and energy. 

Next, we integrate Eq.Q over r and p to obtain rate equations for the particle number. For simplicity, we neglect 
the Hartree-Fock effective potential in the single-particle energy (p,r,t), i.e. 

e,(p,r,t)«e(p,r) = |^ + y. (13) 

We also assume that the temperatures of three internal states are close to each other, so we only retain terms to 
second order in fluctuations of Pi. 

All the terms of the left side of Eq. ([5]) other than the first term do not contribute because those integrand arc odd 
functions. Thus, we have only to integrate Eq. (fT2|) 



= I dv I 77-73^" (P'l-,*)- (14) 



dt J J {2TTny 

Using (fTTj) and (fT^ in and performing integrals explicitly, we obtain the equation of motion for particle number 



as 

dN, 



"•^^i „ r)li \ ^ \ ^ / oa c/3 oa o0 , ca o0 na qP 

ji'j' Q/3 



rN,,Nj, - ^in_^N,Nj \ . (15) 



-(3,,)/2f ' ' [(ft+/3,)/2]' ' 
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FIG. 1: Time evolution of spin populations in a trapped gas of *''Rb atoms. We take the following values for the parameters: 
iV = 1 X 10^, {w^/27r,W;,/27r,w^/27r} = {890,890, 160} Hz and T = 500nK. We set the initial condition as A^i : A^o : A^-i = 
45 : 10 : 45. With this initial condition, the time evolutions of A'^i obtained from Eq. (|17|) and A'^-i obtained from Ea. (|19|) are 
precisely the same, and are indicated by the single solid line. The broken line indicates the time evolution of A^o obtained from 
Eq.lIHJ. 



Here, we defined the coUisional relaxation rate associated with the population transfer as 



2 (3i {muJho) 
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where uj^o = {^x^y^z)^^^ ■ We write down explicit forms of Eq. (|15[) for three components 

'^^l A I AT AT \ n7^ 

dt +/3_i)/2]'' J 
^ = 870 I .A^iiV-i - TVoiVoj , (18) 

^ = 47_. l^N^N, ^'^-^ ^N.nX . (19) 

From Eq. p5p . we find that population dynamics depends not only on the particle numbers but also on the temper- 
atures. In general, the atoms flow into the state with lower temperature. In addition, depending on the ratios of the 
each particle number, spin-spin interaction (32) contributes to increasing mp = particle, or vice versa. 

In general, temperatures vary from time to time. So, we will also derive the equations of motion for the temperatures. 
First, we derive the equation for the energy with neglecting the mean field. The energy of the internal state i is given 

by 

(t) = I dv I (p, r) /, (p, r, t) . (20) 

From the above equation, we can derive the equation for the energy by integrating the kinetic equation ^ multiplied 
by the single-particle energy over real and momentum space. Using the approximations (jlip and (|13p . we obtain 
the expression for the energy in terms of Pi and /i^ 

E, « |- i ^e'^'''-. (21) 
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The energy Ei can also be written in terms of the particle Ni and the temperature parameter (3i as 

With the assumption that the Wigner function obeys the Boltzmann distribution function pip , we obtain the equation 
for the energy as 



dt J J {2TTnf 

Using pTjl and (|12p in (p3)) and performing integrals explicitly, we obtain the equation of motion for the energy as 



'^^i _ o„2 ('^■'^'lo) ^^(qu qI3 Qa aP , Qa qu (^I3 



r)„2 V-'^ho) / qI3 Qa c/3 , QCt qP na q, 

I ji'j' ap 



18[(ft,+/3,,)/2f 9 [ip,+f3,,)/2f 



Eii Eji 



I181(ft + ft)/2f 9|(ft + ft)/2 
Using ([22|) and combing (fT5|) and ([24]) . we obtain the equations of motion for the temperatures 

•^^^ _ r)„2 {'^^hof \^ \^ [ oa qP na qP , na qP qo qP 
^ ' ji'j' ap 

( Pf'Pj' l3^N,.NJ. 5 /3, ,1 A(A> -/3j-0 



[(A,+/3,0/2f 1 6[(ft. +/?,v)/2] 3 +/3^.,)/2]^ 



m+p,)/2f-'[ 6[(ft + /3,)/2] 3p,;+/3^.)/2]^ 



(25) 



We now consider the special situation where the temperatures of three internal states are the same. In this case, 
the temperature stays constant in time. In FiglU we plot the population dynamics in the situation where all the 
temperatures remain the same constant value. We plot the variation of the populations for ^''Rb. where the total 
number of atoms is = 1 x 10'"' and the temperature is T = 500nK, which corresponds to T « 2Tbec- We referred 
several quantities to the experiment of Ref.[13: the trap frequencies are ujz/2n = IGOHz and u;± = 890Hz, and the 
initial situation is Ni : Nq : N-i = 45 : 10 : 45. In FiglU we find that population numbers evolve toward the situation 
where the numbers of particles are the same, where the chemical potentials satisfy Eq.(l9]). This is consistent with the 
behavior of the noncondensate atoms observed in the experiment of Ref. . 

Next, we estimate characteristic time required for the system to reach equilibrium. From Eqs.(|15p and (|16p . we see 
that the characteristic time for the variation of the population is l/7i. The relaxation rate ji is very small because 
7i oc gl, where 52 is the spin-spin interaction (jgol ~ 200 \g2\ for *^Rb) as shown in Eq. fTB)) . A rough estimate gives 
the relaxation time as l/ji x N total — 39s. From the result in Fig[T]for T — 500nK, we determine the characteristic 
relaxation time t as the time at which Ni reach equilibrium value. We find t « 27s. This turns out to be the order 
of the lifetime of trapped atoms. We find that the populations change rapidly in early time when they are far from 
equilibrium. The higher the ratio of each particle number of the internal state, the more rapid the particle numbers 
vary. In the experiment of Ref. [T9l| . the time required for populations to reach equilibrium is t w IDs. This is less than 
half of our result for a thermal gas. Since in the experiment, both the condensate and noncondensate atoms exist, we 
speculate that the interaction between the condensate and noncondensate atoms equilibrates the system much more 
rapidly [1^. We note that the relaxation rate p6)) can be written in terms of the density in the center of the trap 
potential Ui (0): 



In the case of a trapped single-component Bose gas, Ref. (23| showed that in the BEC phase, condensate density 
makes additional contribution to the relaxation rate. Thus, in a trapped gas, the coUisional relaxation time become 
much shorter in the BEC phase because of the high condensate density at the center of the trap. 
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IV. CONCLUSIONS 



111 this paper, we studied the population dynamics of the spin-l Bose atoms above Tbec- We found that the time 
evolution depend on both relative number and relative temperature. Comparing our result with the experiment below 
Tbec J where the condensate is present [Toj . we found that in the case where both the condensate and noncondensate 
atoms exist, equilibrium is reached in less than half the time it takes when condensate atoms are absent. This suggests 
that the condensate atoms help the equilibration of the noncondensate atoms. We also found that the initial time 
evolution of the population number is very rapid nevertheless its relaxation time is very large. 

In future work, we will derive a kinetic theory for spin-l Bose condensed gases below TbeCi where the condensate 
and the noncondensate atoms interact with each other. This theory will consist of coupled equations for the condensate 
and noncondensate. Using this theory, we will study the population dynamics at finite temperatures below Tbec, 
where both condensate and the noncondensate atoms are present. We will also consider coupled collective motion of 
the condensate and noncondensate. 
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